The electronic structures and magnetic properties of perovskite 
ruthenates from constrained orbital hybridization calculations 
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We introduce a method to analyze the effect of hybridization by shifting corresponding atomic 
levels using external potentials. Based on this approach, we study perovskite ruthenates, and 
unambiguously identify that the covalency between the j4-site cation and O ion will modify the 
Ru-O hybridization and change the density of state at Fermi level, consequently affect the magnetic 
properties significantly. We also study the effect of pressure and reveal that hydrostatic pressure has 
a small effect on the Ru-O-Ru bond angle of SrRuOs, while it will decrease the Ru-O length and 
increase the band width significantly. Therefore, the magnetic ordering temperature will decrease 
monotonically with pressure. 

PACS numbers: 71.20.-b, 71.20.Lp, 75.50.Cc 
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Due to the interesting electrical, structural and mag- 
netic properties, perovskite ruthenates have attracted 
continual research attention [ll-[l9|. SrRuOa crystallizes 
in an orthorhombic structure with space group 
Same as SrRuOs, CaRuOs also has an GdFeOs-type per- 
ovskite structure 0, the only observable difference be- 
tween their crystal structures is the Ru-O-Ru bond angle. 
However, their magnetic properties are different com- 
pletely. SrRuOa is a ferromagnetic (FM) metal with a 
Curie temperature Tc=160 K[4|, while there is a debate 
about the ground state of CaRuOaP, [^, [3 [I^- It is 
believed that the leading factor to determine the mag- 
netic properties of these 4rf materials is the structural 
distortion and the associated Ru-O-Ru bond angle (20|. 
Comparing with Sr^+, Ca^+ has smaller radii, therefore 
CaRuOa has bigger structural distortion and smaller Ru- 
O-Ru bond angle, which lifts the band degeneracy and re- 
duces the density of states at the Fermi energy N(E/) (2^ . 
Consequently the Stoner criterion for magnetism may not 
be satisfied for CaRuOa. On the other hand, Ba^+ has 
bigger radii, so perovskite BaRuOa will have larger Ru- 
O-Ru bond angle, comparing with that of SrRuOa. Thus 
it is very natural to expect that BaRuOa should have 
higher T^ than SrRuOs. 



Very recently, using high pressure and high tempera- 
ture techniques, Jin et aZ.[21[, synthesize the cubic per- 
ovskite BaRuOs successfully. The Ru-O-Ru bond angle 
of cubic BaRuOa is 180°, which is indeed larger than that 
of SrRuOa. However, its Tc is only 80 K, which clearly 
indicates that in addition to the bond angle, there is other 
factors which also play important roles in the magnetic 
properties. Jin et aL[2l| propose that the weaker Ba-0 
hybridization will enhance the strength of Ru-O bond, 
consequently reduce the Tc. The importance of ^-site 
cation has also been emphasized by the density func 
tion calculations ^22r,24] . However contrast to Jin 21 



Zayak et a/.[2J] suggest that the magnetization in these 
compounds is anticorrelated with the degree of covalency 



between j4-site ion and O ion. 

In addition to Ru-O-Ru bond angle and ^-site hy- 
bridization, Ru-O bond length is another factor which 
also affect the band structure and the N(Ey), conse- 
quently change the magnetic properties. Hydrostatic 
pressure may vary both the bond length and angle, and 
there are several publications which study the effect of 
pressure p^5t-i28] . All experiments indicate that the Tc of 
SrRuOa will decrease with pressure However, the 

reason why Tc decreases with pressure is still not clear. 
X-ray diffraction[27l suggests the Ru-O-Ru bond angle 
will decrease with pressure, which results in the decrease 
of Tc. On the other hand, neutron powder diffraction 28| 
claims that the Ru-O bond length is sensitive to pressure 
while the bond angle remains almost constant. There- 
fore, to clarify the controversial issue about the influence 
of ^-O covalence and pressure on the magnetic properties 
of perovskite ruthenates is an interesting problem which 
we address in the present work. 

Local spin density approximation (LSDA) is well used 
for perovskite ruthenates (sl [ol. [iol 23, 24, 20|. The main 
aim of our study is to illustrate the effect A-0 hybridiza- 
tion, we thus neglect the effect of Coulomb correlation, 
although there is a debate about whether the electronic 
correlation is important in perovskite ruthenates 0, 
[29l . IsH . We perform the electronic band-structure calcu- 
lations within LSDA using the full potential linearized- 
mufhn-tin-orbital LMTO methodifl|. A mesh of 256 k- 
points in the irreducible Brillouin Zone is used. The self- 
consistent calculations are considered to be converged 
when the difference in the total energy of the crystal do 
not exceed 0.25 meV per formula unit (f.u.) and the dif- 
ference in the total electronic charge do not exceed 10"'^ 
electronic charge at consecutive steps. 

Based on the experimental structure [2l|. we perform 
calculation for the recently synthesized cubic BaRuOa. 
We find that the FM metallic solution is the ground 
state, which agrees with the experiment [21]. The en- 
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ergy difference (Epm-fm) between total energy of its 
PM and FM solution is 25.6 meV/f.u. The magnetic 
moment at Ru site (Mj^„) of BaRuOa is about 0.76 /i^, 
which is also close to the experimental value (about 0.80 
Mb)12i1. Due to the strongly hybridization with the Ru 
4d state, there is considerable moment locating at O site. 
We also perform calculation for SrRuOs and CaRuOs us- 
ing their experimental structure [2]. Consistent with the 
experimental and previous theoretical work, our calcu- 
lation also predict the ground state of SrRuOs is FM 
metal. The obtained Epm-fm is 38.9 meV/f.u., which 
is slightly larger than that of BaRuOa . The magnetic mo- 
ment mainly locates at Ru site (about 0.93 /i^), which 
is also in good agreement with the previous works [io]. 
On the other hand, for CaRuOa, the FM state is unsta- 
ble in energy, which is consistent with other theoretical 
work[20]. 

After reproducing the experimental results, we now try 
to explore the controversy about the effect of A-0 co- 
valence. It is widely accepted that the alkali elements 
at the ^-site of ^BOa perovskite are highly ionic, A- 
site only contribute to the distortion of lattice structure 
and should not affect the Ru-0 hybridization consider- 
ably. However, recently the effect of hybridization be- 
tween A-site cation and O had been emphasized by sev- 
eral groups 2l|-|2J] , although there is debate about its ex- 
act effect. Thus to gain a conclusive comment about the 
effect of A-0 hybridization is important. 

It is well known that the strength of hybridization be- 
tween two orbits is strongly depend on their energy differ- 
ence. We thus introduce a external potential in the Kohn- 
Sham equation to shift the orbital level consequently con- 
trol the hybridization: 

{Hks + \ilma)Vext{ilrna\)i> = E^) 

where, Hks is the Kohn-Sham potential, the basis 
\ilma) is the orbit which we try to shift (i denotes the 
site, n the main quantum number, / the orbital quantum 
number, m the magnetic quantum number and a the 
spin index), and Vext is the magnitude we shift. By this 
constrained calculation, we can directly study how vari- 
ous orbitals and hybridizations between them affect the 
physical properties of solids. Although similar in spirit, 
this constrained orbital-hybridization method is different 
with the constrained LSDA calculation where the exter- 
nal potential is applied to a particular orbital to constrain 
the orbital occupation. 

For perovskite ruthenates, the bands near the Fermi 
energy (E^), which determine the magnetic properties, 
are mainly contributed by Ru Ad and O 2p, on the other 
hand the d band of alkali elements are several eV above 
Ey. Thus when we upshift the empty band of ^-site atom 
by a external potential, the effective hybridization be- 
tween this band and the O 2p band will decrease, there- 
fore we can study the effect of ^-O hybridization directly 
by our constrained orbital-hybridization method. 



Based on the experimental crystal structure^, 21 [, 
we perform constrained orbital-hybridization calculation, 
the numerical results reveal that the position of Sr Ad 
band affects not only the magnetic moment but also the 
magnetization energy Yipm-fm- Upshifting Sr Ad band 
by 9 eV and 20 eV will reduce Mpu and Epm-_fm to 
1.10 20.2 meV; and 0.4 /x^, 7.5 meV, respectively. 
These results indicate the strength of Sr-0 hybridization 
is correlated with the magnetization. Our constrained 
orbital-hybridization calculations also show that upshift- 
ing Ba hd band has a similar effect as in SrRuOa. The 
position of Ca Zd will also affect the band structure of 
CaRuOs, but shifting Ca Zd band cannot make the FM 
state stable. 

To further understand the effect of the ^-O hybridiza- 
tion on magnetic properties, we perform constrained 
orbital-hybridization calculation for PM state and shows 
the density of states (DOS) of BaRuOs in Fig.l. It is 
found that the energy position of Ba bd band has an 
important effect on its electronic band structure. Con- 
ventional LSDA calculation predicts that the Bg band 
of Ru mainly locates at -7.6 to -5.0 eV and -0.2 to 5.1 
eV, while i2g band has two sharp peaks centered at -0.1 
and -5.2 eV, respectively, as shown in Fig. 1(a). Upshift- 
ing Ba 5(i band depresses the hybridization between Ba 
and O, consequently varies the band structure around E/ 
considerably. As shown in Fig. 1(b), decreasing the Ba- 
O hybridization will widen the eg band of Ru slightly, 
while change the i2g band significantly, make the t2g to 
distribute in the whole range from -7.0 to 0.9 eV, and 
decrease N(E/) considerably. It is well accepted that 
^RuOs is an itinerant magnetism, and an large N(Ej) in 
PM calculation is crucial for the magnetization. There- 
fore, according to the Stoner theory, the ^-O hybridiza- 
tion can affect the magnetic properties considerably. We 
also perform constrained orbital-hybridization calcula- 
tion for SrRuOa and CaRuOa. The similar sensitivity of 
electronic structure to energy position of d band of cation 
has also been found. So we can conclude that the ^-O 
hybridization will change the Ru-0 hybridization, conse- 
quently change the magnetization of perovskite ruthen- 
ates. 

To study the effect of pressure on the structure, we 
optimize the structure and relax all independent inter- 
nal atomic coordinates until the corresponding forces 
are less than 1 mRy/a.u. Consistent with the recent 
experiment [i^, we find that the Ru-O-Ru bond angle 
of SrRuOa is not sensitive to the external pressure, while 
the Ru-0 bond length decreases monotonically with pres- 
sure, as shown in Fig. 2. Our numerical results agree with 
the experimental value at low pressures, as shown in the 
inset of Fig. 2, while the experimental bond length at 
high pressure is not available. Based on the optimized 
structure, we perform the band-structure calculation for 
a number of different volumes and fit the curves E(V) 
of the calculated total energies vs volume to the Birch- 
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FIG. 1: (Color online). Partial density of states (DOS) of the 
cubic BaRuOs. Red and green curve denote the t2g and eg 
band, respectively, (a) is the conventional LSDA results; (b) 
is the constrained orbital hybridization calculations results 
with the 5d band of Ba upshifted by 9 eV. 
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FIG. 2: (Color online). Numerical pressure dependence of the 
average Ru-O bond length of SrRuOs (in a unit of A). The 
inset is a comparison between the experimental[23| (red) and 
numerical (black) results at low pressures. 



TABLE I: The obtained equilibrium volume Vo, bulk mod- 
ulus Bo and its pressure derivative Bq of SrRuOa. Previous 
experimental [27I] and the theoretical results from VASP 
and SIESTA package are also listed for comparison. 





Vo(a.u.^) 


Bo (GPa) 




Present LMTO 


1605 


190.5 


4.6 


Experiment [27] 


1633 


192 


5.3 


Other theory (VASP)[31j 




200 


4.6 


Other theory (SIESTA) [SlJ 




219 


4.4 



Murnaghan equation of state. The obtained equilibrium 
volume Vo, bulk modulus at equilibrium Bq, its pressure 
derivative Bq for SrRuOs, together with the available ex- 
perimental and theoretical values are summarized in Ta- 
ble I. Our theoretical equilibrium volume is only about 
1.7% smaller than the experimental values. Such a de- 
viation exists normally in the LSDA calculations. Our 
full potential calculations agree with the experiment very 
well, while the pseudopotential calculation [sH slightly 
overestimates the bulk modulus. 

Perovskite ruthenate belongs to itinerant magnet, 
therefore one cannot estimate the interatomic exchange 
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FIG. 3: (Color online). The numerical pressure dependence 
of the energy difference between PM and FM solution of 
SrRuOa (black), and the experimental,27!| pressure depen- 
dence of the Tc (red). 
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interaction and Tc accurately as in the local magnet [325. 
One important parameter relevant to the Tc is Epm-fm, 
which is plotted in Fig. 3. For comparison, we also plot 
the experimental curve of vs pressure in Fig. 3. It is in- 
teresting to notice that both numerical Epm-fm and ex- 
perimental Tc decrease monotonically and rapidly. The- 
oretical work predicts that the critical pressure (Pc) for 
FM to PM transition is about 21 GPa, after this pres- 
sure FjpM-FM starts to equal or less than zero. The 
experimental Tc at the pressures above 18 GPa is not 
available, but an extrapolation of the experimental data 
at low pressures yields a crude approximation of 23 GPa, 
which is quite close to our numerical results. 

In summary, we introduce a method to identify the 
effect of covalence for a particular material by applying 
orbitally dependent external potential. This approach 
allows us to control the hybridization between the orbits 
directly. Our constrained orbital hybridization calcula- 
tion unambiguously identify that the covalency between 
the ^-site cation and O ion will modify the Ru-0 hy- 
bridization, consequently affect the magnetic properties 
significantly. We also study the effect of pressure and 
find that hydrostatic pressure has a small effect on the 
Ru-O-Ru bond angle of SrRuOa, while it will decrease 
the Ru-0 length, increase the band width and result in 
the monotonic decrease of the magnetic ordering temper- 
ature. 
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